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ABSTRACT 

With the advent of aU-sky Ha surveys it is possible to determine a rehable free-free 
template of the diffuse interstellar medium (Dickinson, Davies & Davis 2003) which 
can be used in conjunction with the synchrotron and dust templates to correct CMB 
observations for diffuse Galactic foregrounds. We have used the COBi^-DMR data 
at 31.5, 53 and 90 GHz and employed cross-correlation techniques to re-evaluate the 
foreground contributions, particularly that due to dust which is known to be partially 
correlated with Ha (and free-free) emission. 

The DMR microwave maps are found to contain, as well as the expected syn- 
chrotron and free-free components, a component tightly correlated to the 140 /im 
dust maps of DIRBE. At 31.5, 53 and 90 GHz this emission is 6.3 ±0.6, 2.4 ±0.4 and 
2.2 ± 0.4 /i,K/(MJy sr^^) at 140 ^m respectively. When corrected for the contribu- 
tion from thermal dust following model 7 of Finkbeiner, Davis & Schlegel (1999), a 
strong anomalous dust-corrclated emission component remains, which is well- fitted by 
a frequency spectrum of the form where (3 ~ 2.5 in the DMR frequency range; 
this is the dominant foreground at 31.5 GHz. The result implies the presence of an 
emission component with a dust-like morphology but a synchrotron-like spectrum. We 
discuss the possible origins of this component and compare it with the recent WMAP 
interpretation. 

The better knowledge of the individual foregrounds provided by the present study 
enables a larger area of the sky (|6| > 15°) to be used to re-appraise the CMB 
quadrupole normalisation, Qrms~PSi and the scalar perturbations spectral index, n. 
We find Qrms-PS= ^^■'^^2 3 with a power-law spectral index of n = 1.2 ± 0.2. These 
values are consistent with previous COBE-DMR analyses and the WMAP 1-year 
analysis. 

Key words: cosmic microwave background - cosmological parameters - radio con- 
tinuum: general - ISM: dust, extinction 



1 INTRODUCTION 

A major goal of cosmology is to determine the fundamen- 
tal cosmological parameters that describe our Universe. Ob- 
servations of anisotropics in the cosmic microwave back- 
ground (CMB) provide unique data in achieving this goal 
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through a determination of their power spectrum. An accu- 
rate measurement of the CMB power spectrum requires a 
good knowledge of the various Galactic foregrounds which 
have their own frequency and power spectra. This is par- 
ticularly important for all-sky, high-sensitivity CMB exper- 
iments such as WMAP (Bennett et al. 2003) and Planck 
(Tauber 2000). For example, the CMB power spectrum ex- 
pected from observations with the Planck satellite will be 
accurate to ~ 1 per cent over the ^-range 100—2000; this cor- 
responds to a temperature accuracy of better than 1 /^K over 
the angular range 2° to 5 arcmin. The Galactic foregrounds 
should therefore be estimated to better than this value, oth- 



2 A. J. Banday et al. 



erwise the CMB data will be degraded. The COBE-HMK 
data (Bennett et al. 1996) covering ^ = 2 — 30 at frequencies 
of 31.5, 53 and 90 GHz are constrained by knowledge of the 
foregrounds; this is also the case for the the WMAP data, 
particularly on the largest angular scales discussed here. The 
challenge is to make use of the largest sky coverage possible 
to reduce the sample variance on such large angular scales 
(iJ 7°); the area available has normally been restricted to 
Galactic latitudes |6| 20° or 30°. In this paper we at- 
tempt to improve on available knowledge of the foregrounds 
and apply it to the CO-B -E-DMR data in order to determine 
whether the fraction of sky available for cosmological assess- 
ment can be extended. The recent WMAP analysis of the 
1st year observations (Bennett et al. 2003a) uses the 23 GHz 
channel as a basis for a Galactic mask which eliminates ~21 
per cent of the sky from detailed CMB analysis (Bennett et 
al. 2003b). 

Previous studies of Galactic foregrounds relevant to 
the CMB at frequencies less than 100 GHz have used the 
408 MHz map of Haslam et al. (1982) as the synchrotron 
template and FIR dust emission (IRAS or DIRBE) as a 
template for the free-free which was thought to be corre- 
lated with dust and Hi (Kogut et al. 1996a,b). It was widely 
established that at 10—50 GHz a dust-correlated foreground 
was present (de Oliveira-Costa et al. 1997; Leitch et al. 1997; 
Kogut 1999) although it was at a level which significantly 
exceeded the free-free emission expected from the partial Ha 
surveys available at the time. This dust-correlated anoma- 
lous emission has been ascribed to spinning dust (Draine 
& Lazarian 1998a, hereafter DL98, 1998b) and appears to 
fit the expected narrow frequency spectrum (de Oliveira- 
Costa et al. 1998, 1999, 2000, 2002), although Mukherjee et 
al. (2001) weaken some of these arguments and claim that 
definitive statements about the origin of the dust-correlated 
component cannot be made based on data at 10 and 15 GHz 
alone. Bennett et al. (2003b) find a dust-correlated compo- 
nent which they interpret as synchrotron emission with a 
flatter spectral index. Lagache (2003), however, argues that 
this emission is intrinsic to the dust itself and arises from 
small, cold grains. This is discussed further in section IFTH 

A significant development now available for assessing 
and quantifying the Galactic foregrounds is the unambigu- 
ous identification of the free-free component independent 
of the quasi-correlated dust. The new Ha surveys by the 
Wisconsin H- Alpha Mapper, WHAM (Reynolds et al. 1998; 
Haffner 1999), and the Southern H- Alpha Sky Survey At- 
las, SHASSA (Gaustad et al. 2001), now make this possible. 
They can be used to produce an all-sky free-free template 
when corrected for Ha absorption by dust and when evalu- 
ated at the electron temperature appropriate to the ionized 
gas in the local region of the Milky Way (Dickinson et al. 
2003). 

In the present paper, using cross-correlation techniques, 
we quantify the synchrotron and free-free emission, leaving a 
clearer estimate of the microwave emission associated with 
Galactic dust. We use the COBE-V)MR data at 31.5, 53 
and 90 GHz (Bennett et al. 1996) to identify an anoma- 
lous dust-correlated component which previously contained 
an unknown amount of free-free emission. Our results are 
strengthened by including correlations with a 19.2 GHz sur- 
vey of the sky (Cottingham 1987; Boughn et al. 1992). 

Our present study leads to a frequency spectral index /3 



for the synchrotron emission; although /3 varies with position 
this global value is adequate for predicting the actual (small) 
contribution at higher radio frequencies {v ^ 30 GHz). We 
also investigate the effects of dust absorption and electron 
temperature in deriving the free-free template. 

Finally, we utilise our improved understanding of the 
three Galactic templates to re-examine the determination 
of the quadrupole normalization Qrms-ps and the scalar 
fiuctuation index n from the COBE-DMK data (Gorski et 
al. 1996) . We extend the area of sky available by using data 
closer to the Galactic plane and examine the stability of 
Qrms-ps and n values as a function of the latitude range 
included in the analyses. 



2 DMR DATA AND THE FOREGROUND 
TEMPLATES 

2.1 Synchrotron template 

The synchrotron emission from the Galaxy originates in 
relativistic cosmic ray (CR) electrons spiralling in Galac- 
tic magnetic fields. These electrons will have a a range of 
spectral energy distributions depending on their age and 
the environment of origin (e.g. supernova explosions or dif- 
fuse shocks in the interstellar medium). The variations in 
CR electron spectral index will translate into variations in 
the frequency spectral index of the synchrotron continuum 
emission. Synchrotron emission, with its temperature spec- 
tral index of /3 = 2.4 to 3.0 (T oc v~^), dominates over 
thermal bremsstrahlung (free-free) emission with (3 = 2.15 
at frequencies below ~ 1 GHz over most regions of the sky. 
Exceptions are the Galactic ridge and the brighter features 
of the local Gould Belt system (Dickinson et al. 2003) 

The only all-sky map at sufficiently low frequency to 
provide a template of Galactic synchrotron emission is the 
408 MHz survey by Haslam et al. (1982). This has a reso- 
lution of 51 arcmin and has a brightness temperature scale 
"absolutely calibrated" with the 404 MHz 8°.5 x 6°.5 survey 
of Pauliny-Toth & Shakeshaft (1962). 

In the present investigation we will derive (for the sky 
area analysed) a mean spectral index between 408 MHz and 
19.2 GHz from the Cottingham (1987) survey and 31.5, 53 
and 90 GHz from the COBE-DMR survey. It is important 
to remember that the spectral index of Galactic synchrotron 
emission varies with frequency and position on the sky. The 
analysis of low- frequency surveys by Lawson et al. (1987) 
indicates that the brighter regions away from the Galactic 
plane have typical spectral indices of /3 = 2.55 and 2.8 at 100 
and 800 MHz respectively. Furthermore, the value of /3 eval- 
uated between any pair of frequencies separated by a factor 
of 2 — 3 varies typically by A/3 — ±0.15 over the sky (Reich 
& Reich 1988). Bennett et al. (2003b) propose a model us- 
ing WMAP data which includes a flat spectrum synchrotron 
component (/3 2.5) in regions of star formation near the 
Galactic plane and a steep spectrum component (/3 ~ 3.0) 
from the Galactic halo. 

At frequencies higher than 2 — 3 GHz no reliable zero 
levels can be measured for large-area surveys at intermedi- 
ate and high latitudes where CMB measurements are nor- 
mally conducted and as a consequence, large area spectral 
index maps cannot be made. However, the spectral index 
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of smaller angular scale (~ 5° and less) structures can be 
estimated from cross-correlation analyses of maps at differ- 
ent frequencies. For example, Hancock et al. (1997) found a 
spectral index of /3 = 3.0 at 10 GHz from bremsstrahlung 
data at 10,15 and 33 GHz. It is interesting to note that the 
steepening of the brightness temperature spectral index by 
0.5 at higher frequencies expected from synchrotron losses 
would lead to /3 = 3.1 at higher frequencies; this is the value 
at these frequencies expected from the local CR electron 
spectrum (Platania et al. 1998). 

2.2 Hq free-free template 

For the free-free template, we use the all-sky Ha template 
described in Dickinson et al. (2003). The all-sky map is a 
composite of Ha data from the WHAM survey of the north- 
ern sky (Reynolds et al. 1998; Haffner 1999) and SHASSA 
survey in the southern sky (Gaustad et al. 2001). The 
SHASSA data are first smoothed to the WHAM resolu- 
tion of 1° and are preferred at declinations 5 < —15°. It 
should be noted that the all-sky template is a first attempt 
at producing a high-sensitivity all-sky Ha map. SHASSA 
data have significantly lower sensitivity than data from the 
WHAM survey; the noise levels are estimated as 0.5 and 0.05 
Rayleigh {R) respectively, where 17?= 10®/47r photons s"^ 
cm"^ sr"^ = 2.41 x lO'^ erg s"^ cm"^ sr"^ = 2.0 cm"*'pc 
for Te = 7000 K gas. Baseline uncertainties exist in both 
data-sets, particularly SHASSA data which results in loss of 
information on angular scales ^ 10° for data at 5 ^ —30° 
due to the 13° field-of-view of each pointing. However, the 
SHASSA background levels are fitted to WHAM data where 
WHAM data exist [5 > —30°) with an assumed cosecant 
law for declinations further south. The different observing 
strategies result in different beam shapes which have not 
been accounted for. Finally, large Ha signals (emission and 
absorption) from bright stars have been filtered using differ- 
ent techniques; the WHAM team use data from nearby unaf- 
fected pixels, while the SHASSA team use a median filter to 
reduce star residuals. However, all these possible systematics 
are significantly reduced when the data, at a resolution of 
~ 1°, are convolved with the DMR beam at ~ 7° resolution. 
Indeed, we have explicitly verified this by comparing the 
power spectra of the high-latitude Ha emission after vary- 
ing the latitude at which the surveys are matched and ad- 
justing baselines. Recently, Finkbeiner (2003) has produced 
an all-sky Ha map at 6 arcmin resolution which includes 
higher resolution data from the Virginia- Tech Sky Survey 
(Dennison et al. 1998). The differences in the power spec- 
trum of our Ha map with that given by Finkbeiner (2003) 
are less than 10 per cent for ^ < 15 and less than 30 per cent 
everywhere. 

The optical Ha line (A = 656.3 nm) is well known to be 
a good tracer of the free-free continuum emission at radio 
wavelengths since they are both emitted by the same ion- 
ized gas. Since they are both proportional to the Emission 
Measure {EM = J nidi) there is a well-defined relationship 
between the Ha line intensity /hq and the radio contin- 
uum brightness temperature (Osterbrock 1989; Valls- 
Gabaud 1998; Dickinson et al. 2003). There is a strong de- 
pendence on frequency since free- free emission scales as v^'^^ 
at GHz frequencies. There is also a modest dependence on 
the electron temperature (oc T^'''). In the Galaxy, Te can 



range between 3000 K and 15000 K (Shaver et al. 1983); the 
majority of Hll regions have Te between 5000 and 10000 K 
(Reynolds 1985). Shaver et al. (1983) show that the mean 
electron temperature of Hll regions at the galactocentric ra- 
dius of the Sun {Ro = 8.5 kpc) is 7200 K with an r.m.s. 
spread of 1200 K. In the region of the Galaxy covered in 
the present study where |&| > 15° the ionized hydrogen will 
be typically within ~ 1 kpc of the Sun and accordingly we 
adopt a value of Te = 7000 K. With this value of Te, the 
conversions from Ha to free-free continuum for frequencies 
31.5, 53 and 90 GHz are 5.25, 1.71 and 0.54 yK R'^ respec- 
tively. In section^Jwe will use this prior information to allow 
a fixed free-free component to be subtracted giving a clearer 
separation of the foregrounds. We will take a variation of 
±2000 K from the local value of 7000 K as in indication of 
the upper and lower bounds allowed. Other caveats become 
important at the few per cent level: i) the helium abundance 
is usually taken to be 8 per cent with a spread of a few per 
cent ii) the emitting medium is assumed to be optically thick 
(case B) to the Lyman continuum photons. 

A major uncertainty when using the Ha template is 
the absorption of Ha by foreground dust. The absorption 
can be estimated using the 100 nm maps given by Schlegel, 
Finkbeiner & Davis (1998, hereafter SFD). The absorp- 
tion at Ha A(Ha) in magnitudes using their temperature- 
corrected column density map (D^) in units of MJy sr~^ 
is then A{Ha) = (0.0462 ± 0.0035)L»^/d where fa is the 
fraction of dust tn front of the Ha actually causing absorp- 
tion in the line of sight. The corrected Ha intensity is then 
■^H " = X ioO''xO-4>!/<jxo.0462 Although there are uncer- 
tainties in this relationship due to different reddening values 
expected from different dust populations, the largest uncer- 
tainty is in fd, the fraction of dust lying in front of the 
Ha-emitting region. Dickinson et al. (2003) show that for a 
range of longitudes {I = 30° — 60°), and for |fe| = 5° — 15°, 
fd ~ 0.3, while for local high latitude regions such as Orion 
and the Gum Nebula, there is little or no absorption by 
dust {fd ~ 0). For much of the high latitude sky, is ^ 5 
MJy sr~^ and therefore the absorption is likely to be small 
(< 0.2 mag). At lower Galactic latitudes |fe| < 10°, the dust 
absorption strongly increases to 1 mag and above. For the 
majority of the Galactic plane {\b\ < 5°), the absorption 
is too great to make any reliable estimate of the true Ha 
intensity. 

Finally, the scattering of Ha light from dust at high 
Galactic latitudes could, in principle, cause a degeneracy in 
the analysis since the scattered Ha light would be directly 
correlated with the dust. However, Wood & Reynolds (1999) 
predict that this accounts for less than 20 per cent of high 
latitude Ha emission. Since the anomalous component is 
larger than the free-free component at frequencies higher 
than ~ 20 GHz (see sections 0| and |^ , this systematic will 
not significantly affect our results within the errors. 



2.3 Dust templates 

As templates for the thermal and anomalous dust-correlated 
emission, we utilise the CO-B -E-DIRBE data repixelised into 
the HEALPix pixelisation scheme. This was achieved us- 
ing the publicly available Calibrated Individual Observa- 
tions (CIO) files together with the DIRBE Sky and Zodi 



4 A. J. Banday et al. 



Atlas (DSZA, both products being described in Hauser et 
al. 1998) , allowing the construction of DIRBE full sky maps 
corrected for the zodiacal emission according to the model 
of Kelsall et al. (1998). In particular, we use the DIRBE 
140 /im data at 0°.7 resolution: although the 100 /^m data 
has the best sensitivity, there remain more noticeable arti- 
facts present in the sky map after zodiacal correction. While 
this is not likely to present a problem after the data has been 
convolved with the DMR beam profile (~ 7° FWHM), our 
choice allows maximal consistency with previous analyses 
by the DMR team (Banday et al. 1996; Gorski et al. 1996; 
Kogut et al. 1996a, 1996b). Furthermore, the SFD version 
of the 100 /im sky map has been used to correct for dust 
absorption in the Hq template, thus our choice allows some 
additional independence in the dust correlation results. 

DL98 have suggested that the anomalous dust- 
correlated signal which we re-investigate in this work should 
be significantly correlated with 12 /im emission, so that the 
corresponding DIRBE sky map may prove a better tracer 
for this foreground component. However, the 12 yLxn map 
contains residual zodiacal artifacts in the region of the eclip- 
tic plane which remain significant even after smoothing to 
DMR resolution. The impact of this could be minimised by 
excluding a region close to the ecliptic plane, but we prefer 
to make use of the greatest sky coverage possible. Further, 
de Oliveira-Costa et al. (2002) have tested this hypothesis, 
but it is clear from their fig. 4 that even after considerable 
processing of the 12 /xm DIRBE data for the removal of both 
zodiacal emission and point sources, the longer wavelength 
data at 100, 140 and 240 /xm remain equally viable for cor- 
relation studies. 

Finally, in order to attempt to isolate the anomalous 
dust-correlated component signal, it is necessary to identify 
and subtract the thermal dust emission. We have employed 
the predictions of Finkbeiner, Davis & Schlegel (1999, here- 
after FDS) for the thermal dust emission at DMR wave- 
lengths. EDS consider 8 models based on the 100 and 240 /xm 
maps from SFD, and test them against the COBB-FIRAS 
data. The inconsistency of extrapolations of the SFD maps 
using single component power-law emissivity functions with 
the FIRAS data over the range 200 - 2100 GHz motivated 
the development of generalised two-component models. We 
have considered FDS models 7 and 8 as plausible tracers 
of the thermal dust component emission. Masi et al. (2001) 
have demonstrated that model 8 extrapolated to 410 GHz 
is well correlated with the BOOMERanG data at that fre- 
quency, and a recent analysis of the MAXIMA data (Jaffe et 
al. 2003) also see consistency with these extrapolated dust 
models. It should be remembered that the FDS models are 
formally correct only on scales comparable to or larger than 
the FIRAS beam, hence the FDS predictions of model 7/8 
should be well-matched to our purposes. 

2.4 COBB-DMR data 

In this analysis, we utilise the COBE-DMK 4- year sky maps 
(Bennett et al. 1996) at 31.5, 53 and 90 GHz in Galactic 
coordinates and the HEALPix^ pixelisation scheme with a 
resolution parameter of A^sidc = 32, corresponding to 12 288 

^ http://www.eso.org/science/healpix/ 



pixels on the sky. Inverse-noise-variance weighted combina- 
tions of the A and B channels are formed at each of the three 
DMR frequencies, the weights being (0.618, 0.382), (0.579, 
0.421) and (0.382, 0.618) at 31.5, 53 and 90 GHz respec- 
tively. The dominant feature visible in the sky maps (given 
that the dipole anisotropy of amplitude ~ 3 mK is largely 
removed in the pre-map-making data processing) is the en- 
hanced emission along the Galactic plane. This emission is 
indeed a combination of the synchrotron, free- free and dust 
emission that we are interested in, however, the distribution 
and spectral behaviour at the lowest latitudes {\b\ < 10°) 
is sufficiently complex that we do not attempt to fit or 
model it, considering only the mid- and high-latitude Galac- 
tic emission. Contrary to previous analyses, we attempt to 
fit our templates to somewhat lower latitudes, the results of 
which are discussed in section 2] 

2.5 19.2 GHz data 

In addition to the COBE-T)MK data, we also consider cor- 
relations of the foreground templates with a 19.2 GHz sur- 
vey of the sky (Cottingham 1987; Boughn et al. 1992). Such 
an investigation is clearly of interest here since 19.2 GHz is 
close to the frequency at which the DL98 models for spin- 
ning dust show a maximum in spectral behaviour. Even if 
the anomalous dust-correlated component is not the puta- 
tive rotational dust emission, the coupling amplitude of the 
dust template at 19.2 GHz should enable the spectrum of 
this foreground to be established, in particular since the fre- 
quency of observation is well-removed from that where a 
significant thermal dust contribution is present. The data 
are available in the COBE QuadCube format (White & 
Stemwedel, 1992) with 24,576 pixels on the full sky of size 
1°.3 X 1°.3, and contains observations of angular resolution 
3° FWHM. A small region close to the South Ecliptic pole 
was not observed. The typical noise level is quoted as 1.5 
mK per resolution element, but there is both variation in 
the number of observations per pixel and pixel-pixel covari- 
ance. In contrast to the analysis of de Oliveira-Costa et al. 
(1998) we take account of the former, but the information is 
not available for the latter. Some attempt to provide addi- 
tional constraints on the noise uncertainty is performed by 
scaling the noise level by a factor in the range 0.9 - 1.1, and 
including this parameter in our later statistical analysis. 



3 ANALYSIS METHOD 

Each of the three DMR maps, the 19.2 GHz data, and the 
foreground sky map templates (after convolution to the cor- 
rect angular resolution) are decomposed into harmonic co- 
efficients using a basis of orthonormal functions explicitly 
computed on the incompletely sampled sky (Gorski 1994), 
a procedure which is directly analogous to the expansion 
of a full-sky map in terms of spherical harmonic coefficients. 
This technique has the advantage of allowing exact exclusion 
of the monopole and dipole components from the analysis 
(which are unimportant for the investigation of cosmological 
anisotropy), and makes full use of all of the available spatial 
(phase) information. 

We have adopted several definitions for the area of sky 
coverage, in particular the standard Galactic cut defined as 
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the region with Galactic latitude, above 20° with addi- 
tional regions at Ophiuchus (up to b « 30° at Z ~ 0°) and 
Orion (down to 6 « -35° at / ^ 180°) based on the DIRBE 
140 map (Banday et al. 1997, but recomputed explicitly 
for the HEALPix scheme), together with more aggressive 
and conservative cuts of |6| greater than 15° and 30° respec- 
tively. For each sky coverage a new basis of orthonormal 
functions is computed. This procedure must also be per- 
formed separately for the 19.2 GHz data due to its different 
pixelisation scheme and pixel resolution (and indeed for the 
foreground templates matched to the same pixelisation and 
angular resolution). 

The measured map coefficients in the new orthogonal 
basis at a given frequency, , can then be written in vector 
form as Cv = cqmb + cn -f oixcx where cqmb, cn and 
cx are the coefficients describing the true CMB distribution, 
the noise and one foreground template map respectively, ax 
is now a coupling constant (with units /iK X~^) to be de- 
termined by minimizing 

= (ci^ ^ cexcx)'^ M ^ ~ axcx) (1) 

where M is the covariance matrix describing the correla- 
tion between different Fourier modes on the cut-sky and is 
dependent on assumed CMB model parameters and the in- 
strument noise. Full details of its computation can be found 
in Gorski (1994) and Gorski et al. (1996). 

By minimising the with respect to ax we find that 
the coupling constant ax has an exact solution 

ax — Cx M c,/cxM cx (2) 
and a Gaussian error given by 

a'(ax) = (cxM"'cx)-' (3) 

The contribution of component X to the observed sky 
is the map axTx (where Tx is in the natural units of the 
foreground template map so that this scaled quantity is in 
HK). 

In this analysis, however, we treat multiple foreground 
templates and all three DMR channels simultaneously. This 
forces the CMB distribution (in thermodynamic units) to be 
invariant between the three frequencies, and suppresses the 
sensitivity of the method to noise features in the individual 
DMR maps which are clearly not common to all three fre- 
quencies, thus allowing an improved determination of the ex- 
tent of chance cosmic or noise alignments with the template 
maps. We have tested the impact of such chance alignments 
to the method by repeating the fits after rotation of the fore- 
ground templates. We only find significant correlations when 
the templates are correctly aligned. Moreover, the r.m.s val- 
ues of the coupling coefficients derived from the rotated sky 
maps are consistent with the errors determined directly in 
the analysis. The full details of this method are contained 
in Gorski et al. (1996). The 19.2 GHz data are treated in 
a similar fashion, but we do not extend the method to in- 
clude it in a simultaneous fit with the DMR data to the 
foreground templates because of complications introduced 
by its different angular resolution and noise properties. 

Having solved for the values of the coupling coefficients 
and their Gaussian errors, we are then in a position to com- 
pute the likelihood of the foreground corrected data as a 
function of a given set of input cosmological parameters as 



contained in the matrix M of equation Indeed, this also 
allows us to assess the sensitivity of the coupling constant 
fits to the assumed cosmological power spectrum. We con- 
sider a full grid of power law models specified by the rms 
quadrupole normalisation amplitude Qrms~ps and a spec- 
tral index n as in Gorski et al. (1996). The cosmological 
implications of the foreground corrections as reappraised in 
this work are discussed in section |S| 

In order to quantify the uncertainty in the Hq template 
resulting from dust absorption, we generate corrected Ha 
templates for three values of fd - 0.0, 0.5 and 1.0 - which are 
then fitted successively to the data over the full Qrma-ps, 
n grid of cosmological models, and also test templates cor- 
rected assuming values of fd running from to 1 in steps of 
0.01 for a scale-invariant spectrum with Qrms~ps ~18 /iK. 
Since the correction is non- linear in fd (see section 12.21 . 
modifications to the template are made before convolution 
with the appropriate beam response, at the 1° effective an- 
gular resolution of the Ha template (see section for full 
discussion) . 

We also analyse the DMR and 19.2 GHz data after cor- 
rection for free-free emission using the Ha templates scaled 
by the factors given in Table 4 of Dickinson et al. (2003). 
This enables us to test the impact of cross-talk between the 
templates, which is likely to be present at some level due to 
the large-scale cosecant variation with latitude of the Galac- 
tic emission present in each template (see section . 

Finally, in an attempt to unambiguously identify the 
anomalous dust-correlated emission, we use the FDS mod- 
els of thermal dust emission and subtract them from the 
observed data before correlation studies are performed (see 
section 



4 RESULTS - FOREGROUNDS 

The results of our correlation study for the three DMR fre- 
quencies are summarised in Figs n8')-(c). This summary 
gives information for six parameters as foUows:- 

1. Frequency. Data are grouped in each of the three 
figures (a)-(c) as three bands at each of the DMR frequencies 
31.5, 53 and 90 GHz. 

2. Foreground components. The three panels in each 
figure give the results for dust, free-free and synchrotron 
plotted from top to bottom. 

3. Galactic cut. The results are shown for three Galactic 
cuts, namely |6| > 15° (squares), the standard Galactic cut 
of Banday et al. (1997) (circles) and |&| > 30° (diamonds). 

4. Ha extinction. The effect of the three levels of Ha 
extinction, fd = 0,0.5 and 1.0, are shown in Figs (a), (b) 
and (c) respectively. 

5. Forced/unforced fits to the Ha data. Open symbols 
are fits that are made to the dust and synchrotron templates 
after correcting the DMR data for the free- free contribution 
predicted by the Ha template; filled symbols correspond to 
fits made to the three templates simultaneously in which no 
constraints have been placed on the Ha template correlation. 

6. Electron temperature. The free-free predictions from 
the Ha template depend upon the electron temperature; the 
dashed boxes in the middle panels show the range of pre- 
dicted free-free emission for temperatures 5000 to 9000 K 
(open symbols are for 7000 K). 




31.5 53 90 31.5 53 90 31.5 53 90 

^ 0.0 Frequency [GHz] ^ 0.5 Frequency [GHz] ^ 1.0 Frequency [GHz] 

Figure 1. Derived coupling coefficients between the DMR data at 31.5, 53 and 90 GHz and three templates for foreground emission: dust is traced by the DIRBE 140 fim sky map 
(top panels), free-free emission by Ha data (middle panels), and synchrotron emission by a sky map at 408 MHz (bottom panels). Coefficients are computed for three Galactic cuts: 
|6| > 15° (squares), the standard Galactic cut (circles) from Banday et al. (1997) and |6| > 30° (diamonds). The filled symbols correspond to fits made between the DMR data and the 
three templates simultaneously, while open symbols are fits made to the dust and synchrotron tracer maps after correcting the DMR data by the free-free emission predicted by the 
Ha template according to the scalings in Dickinson et al. (2003). In the upper and lower panels, the closed boxes show results consistent with previous correlation studies using only 
the DIRBE and Haslam templates and the standard Galactic cut. In the middle panel the dashed boxes show the coupling coefficient amplitudes predicted by Dickinson et al. (2003) 
for the Ha -correlated free-free emission assuming values for the electron temperature between 5000 K and 9000 K (the lower and upper bounds respectively). Figs a - c are results for 
three assumed values (0,0.5 and 1.0) of the dust absorption parameter (see section 2.2). 



Reappraising the COBE-DMR data 7 



The results obtained from an analysis using only the 
DIRBE and synchrotron templates, consistent with previous 
work by the COBS-DMR team (Kogut et al. 1996b), are 
shown with their errors by full-line boxes in the upper and 
lower panels of each figure^. This allows a direct assessment 
of the impact of the inclusion of a free-free template on the 
analysis by comparison with results in which such a template 
has not been utilised. 

The results presented here have been derived assuming 
a scale-invariant spectrum (n = 1) for the CMB contribu- 
tion with an r.m.s quadrupole normalisation of Qrms-ps 
— 18 fiK. However, as part of our studies of the impact of 
the foreground correction for cosmological interpretations of 
the DMR data, we have determined the coupling parameters 
over a grid in Qrms~ps and n. We compute an average of 
the a coefficients weighted by the likelihood function over 
the grid and find that the marginalised coupling coefficients 
are essentially identical to those shown here, and that no 
conclusions derived from later analysis are compromised by 
adopting these values as canonical. As described above, we 
have adopted three Galactic cuts, and also computed corre- 
lations in which the coupling to the Ha template has been 
left unconstrained, or where we have fixed and subtracted 
the predicted free- free contribution according to the spectral 
behaviour predicted in Dickinson et al. (2003), determined 
for an electron temperature of 7000 K, the mean value in 
the solar neighbourhood (see section TITIl . 

The following sections provide a detailed discussion of 
this data set. 



4.1 DMR and Dust 

In order to determine accurate values of the coupling am- 
phtudes between the DMR data and the three foreground 
templates it is important that the templates are uncorre- 
lated or largely uncorrelated, or that we know the precise 
emission expected from one of the partially correlated tem- 
plates. The extent of the partial correlation between the 
dust and the free-free (Ha) emission can be elucidated in 
the present study by first removing the free-free emission 
derived from the Ha template (section 12.211 before deriving 
the coupling coefficient for DIRBE dust. The computation 
was also made without this Ha constraint; the results were 
not significantly different. It is also necessary to correct the 
Ha template for absorption by intervening dust by an ab- 
sorption factor fd; the analysis of Dickinson et al. (2003) 
indicates that this is between and 0.3. At high Galactic 
latitudes this has little effect since the total absorption is 
small (A(Ha< 0.1). The plots in Fig. [T] show the effect of 
assuming different values of fd- In fact, the stability of the 
dust couplings derived when the free-free coupling is uncon- 
strained implies that there is minimal cross-talk between the 
components. The results when the free-free coupling is fixed 
are more strongly infiuenced by the relation between Te and 
fd as discussed in more detail in section 

^ Note that the boxes in the top and bottom panels of Fig.lTlshow 
slightly different results to the previous analysis of Kogut et al. 
(1996b). However, this is not surprising since they used an earlier 
release of the DIRBE 140 fim map and the data were analysed in 
the QuadCube format. 



The estimates of the coupling amplitudes of DMR data 
against the DIRBE 140 fim template at the the three fre- 
quencies of DMR are given in the top panels of Fig. The 
best estimates are for fd~0 (see discussion in section [4.31 . 
for the largest sky area covered (|6| jj 15°) and for a fixed 
Ha contribution. The values are 6.3 ±0.6, 2.4 ±0.4, 2.2 ±0.4 
/iK/(MJy sr"^) at 31.5, 53 and 90 GHz respectively. We 
compare these values with those of Kogut et al. (1996b) us- 
ing only the DIRBE and the Haslam templates of 6.4 ± 1.5, 
2.7 ± 1.1 and 2.8 ± 1.0 /iK/(MJy sr"^). There is a smaU re- 
duction in coupling amplitude when the free-free (Ha) tem- 
plate is taken into account. Significantly, there is a factor of 
2 — 3 reduction in the r.m.s. error when the free-free emis- 
sion is included in the analysis. This reduction is partly due 
to the larger area of sky covered but is mainly due to the 
substantial scatter in the correlation between dust and Ha 
emission (Dickinson et al. 2003) and the fact that they are 
individually cross-correlated with the DMR data. 

We will now consider the r.m.s. brightness temperature 
contribution expected for each foreground component as cal- 
culated from the correlation coefficients in Fig.0 The entries 
in Table are the product of the correlation coefficient and 
the r.m.s. level on each of the template maps (after subtrac- 
tion of the best-fitting monopole and dipole) for the relevant 
range of Galactic latitudes and for each of the three DMR 
frequencies. In Tabled we have converted the conventional 
(Rayleigh- Jeans) brightness temperature units into thermo- 
dynamic temperature units (relative to a fixed CMB tem- 
perature, Tcmb) by multiplying by the Planck correction 
factor {e^ — 1)'^ /{x'^e^) where x = hv /kTcms', the factors are 
~ 1.03, 1.07 and 1.23 at 31.5,53 and 90 GHz respectively. At 
the preferred value of fd~0, the rms dust emission at 31.5 
GHz for the |6| > 15° cut is ~ 4 times that for \b\ > 30° re- 
flecting the strong increase towards the Galactic plane. The 
table shows that the dust emission is ~ 1.5 times stronger 
than the free-free, the next brightest component. We note 
that the r.m.s. error in the dust foreground estimate is sim- 
ilar (-^ 6 fiK) for each Galactic cut; this suggests that the 
precision with which the CMB can be corrected for dust 
emission is not degraded on extending the field of view to 
|6| ~ 15°. 



4.2 DMR and Synchrotron emission 

Although synchrotron emission is weak at DMR frequen- 
cies, some significant conclusions can be drawn from the 
cross-correlation between the DMR data and the 408 MHz 
synchrotron template. The most compelling data are shown 
in the bottom panel of Fig. 0^a) which is for fd = 0.0, and 
for the cross- correlation obtained with the free-free emis- 
sion predicted from the Ha template subtracted, as shown 
by the open squares. The couphng constants ^^Synchrotron 
for the largest sky coverage (|6| > 15°), are 0.91 ± 0.83, 
-0.26 ± 0.56 and -0.35 ± 0.54 ^iK K'^ at 31.5, 53 and 90 
GHz respectively. The coupling constant can be converted 
to limits on the brightness temperature spectral index (5 
between 0.408 GHz and the observing frequency, v, using 

Synchrotron 

^ V 0.408 

Table |5]gives the values of [3 for the three sky cuts used 
in the present study. The most significant results are for 
the |b| > 15° cut which includes the brightest backgrounds; 
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Table 1. The predicted foreground r.m.s signal in the DMR sky maps at ^ 7°resolution due to dust, free-free and synchrotron emission, 
determined by scaling the corresponding emission templates by the best-fit coupling coefficients. Quoted errors are 68% confidence level, 
all values are in /xK and have been converted to thermodynamic temperature relative to a fixed CMB temperature, Tcmb = 2.73 K . We 
do not distinguish between the thermal and anomalous dust-correlated contributions. Note the lack of sensitivity to for the free-free 
contribution. 



Component Coverage DMR Map (GHz) 









31.5 






53 






90 








fd = 0.0 


0.5 


1.0 


fd =0.0 


0.5 


1.0 


fd =0.0 


0.5 


1.0 




b > 15° 


51. 8± 6.0 


41. 6± 6.6 


42. 2± 6.6 


18.6± 3.7 


14.8± 4.1 


15.4± 4.0 


19.6± 4.3 


18.6± 4.8 


19. 6± 4.7 


Dust 


standard 


19.9± 5.9 


17.9± 6.1 


17.3± 6.3 


lO.Oit 4.1 


8.9± 4.3 


8.7± 4.3 


ll.Oit 4.5 


9.9± 4.7 


9.6± 4.8 




b > 30° 


13.6± 7.1 


10.6± 7.4 


9.9± 7.6 


14.0± 4.8 


11. 8± 4.9 


11. 7± 5.1 


10.8± 5.3 


8.8± 5.5 


8.3± 5.7 




b > 15° 


29.1± 4.4 


35.8± 4.8 


35. 4± 4.8 


12. 9± 2.6 


14.2± 2.8 


13.8± 2.7 


5.3± 3.0 


4.1± 3.3 


3.4± 3.3 


Free-free 


standard 


lO.Oit 4.2 


11. 2± 4.2 


11. 5± 4.3 


2.9± 2.5 


3.5± 2.5 


3.4± 2.5 


3.6± 2.9 


4.1± 2.9 


4.2± 3.0 




b > 30° 


10.8± 5.6 


13.6± 5.4 


14.1± 5.7 


-3.9± 3.5 


-0.9± 3.1 


-0.8± 3.2 


0.3± 4.0 


2.7± 3.8 


3.2± 3.9 




b > 15° 


5.3± 5.4 


8.5± 5.4 


8.5± 5.5 


-1.6± 3.8 


-0.5± 3.8 


-0.6± 3.8 


-2.7± 4.2 


-2.5± 4.2 


-2.6± 4.2 


Synchrotron 


standard 


4.7± 5.4 


5.0± 5.4 


5.2± 5.4 


0.2± 3.7 


0.2± 3.7 


0.3± 3.7 


-3.2± 4.1 


-3.1± 4.1 


-3.0± 4.1 




b > 30° 


8.0± 5.9 


8.6± 5.9 


8.7± 5.9 


-4.5± 3.9 


-4.5± 3.9 


-4.5± 3.9 


-2.7± 4.4 


-2.6± 4.4 


-2.5± 4.4 



Table 2. The derived synchrotron spectral index between 
408 MHz and the quoted frequency, assuming the power-law re- 
lation a^yn-chrotrori ^ (^_^~)-'' Quoted errors are 68% con- 
fidence level, upper limits are at the 95% confidence level. The 
numbers are determined from the coupling amplitudes resulting 
from an analysis adopting a fixed free-free contribution. 



Frequency (GHz) 




Sky Coverage 






|6| > 15° 


Standard 


|fe| > 30° 


19.2 


3.08t«:J° 


o nn+O.lO 
•^•"^-0.07 


>3.08 


31.5 


>3.07 


>3.01 


>2.93 


53.0 


>2.83 


>2.78 


>2.72 


90.0 


>2.58 


>2.54 


>2.49 



~ 95% (2cr) upper limits are given for /?. The 19.2 GHz 
result will be discussed in section ITU 

The spectral index calculated between 0.408 GHz and 
20 to 30 GHz where synchrotron signals can be detected 
is /3 = 3.1. /3 is known to increase with frequency between 
~ 0.1 and 1 GHz (see section 12. H with typical values of 
2.7 at 1 GHz. Clearly /3 steepens further by 20 - 30 GHz; 
the "local" value of the index will be greater than 3.1 at 
DMR frequencies if the synchrotron-loss turnover (A/3 = 
0.5) occurs at frequencies higher than a few GHz. Spectral 
indices of ~ 3 were considered by Banday & Wolfendale 
(1991). Direct measurements at frequencies of 5, 10 and 15 
GHz in the Tenerife and Jodrell Bank experiments (Jones 
et al. 2001, Hancock et al. 1997) indicate spectral indices 
of ~ 3 between 0.408 GHz and these frequencies. We note 
that the values derived in Table |21 are the average spectral 
index over a substantial part of the sky; there is presumably 
a scatter of spectral index at high frequencies corresponding 
to that found at low frequencies. 

Table shows that the r.m.s. synchrotron emission in 
the DMR data is essentially the same for each of the three 
latitude cuts in contrast to the dust and free-free emis- 
sion which increases significantly towards the Galactic plane. 
This is because in the 408 MHz synchrotron template map. 



the North Polar Spur (NFS) dominates the emission in each 
latitude cut. The NFS is the major contribution to the syn- 
chrotron r.m.s. temperature at these frequencies. This will 
be discussed further in section 



4.3 DMR and Free-free 

From Fig. Q the Ha template appears to be well correlated 
against the DMR data at all frequencies for the 15°Galactic 
cut, and at better than the 2a level at 31.5 GHz for each of 
the cuts. The derived correlations with the Ha template are 
completely consistent with that expected for free-free emis- 
sion due to thermal electrons with temperatures in the range 
5000 — 9000 K. For the larger values of fd, the picture is less 
consistent, and it appears that some of the Ha correlation is 
at the expense of the dust, particularly for the case of fixed 
scaling. For a value of fd of order unity, the derived scaling 
for a 15°cut is inconsistent with this temperature range. 

In order to study the dependence on the parameter fd 
more extensively, we have run a further set of fits to the 
data in which we keep Qrma-ps and n fixed at values of 
18 and 1 respectively, and vary fd over the range - 
1 in steps of 0.01. In Fig. HI the results for the 31.5 GHz 
channel over the three Galactic cuts of choice in our analy- 
sis are shown. As before, it is clear that the derived scaling 
amplitude decreases as fd increases. This has a simple phys- 
ical interpretation, namely that the amount of foreground 
emission correlated with the Ha template as admitted by 
the 31.5 GHz sky map is fixed (see Table 0, and since the 
predominant effect of the correction for extinction is simply 
to enhance the Ha signal with very little change in detailed 
structure at high latitudes and over the range of scale probed 
by CO-B -E-DMR, the coupling amplitude must fall with in- 
creasing fd as observed. This observation then allows us to 
restrict our investigations to low values of fd at high lati- 
tudes, consistent with the results of Dickinson et al. (2003). 

In fact, we can take the analysis yet further. By noting 
that the scaling amplitude is effectively related to the value 
of the electron temperature in the interstellar medium (Te), 
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Figure 2. The coupling coefficient determined between the DMR 
data at 31.5 GHz and the Ha tracer of free- free emission as a 
function of the dust absorption parameter f^. The grey bands 
correspond to the 95% confidence ranges for the derived cou- 
pling coefficient amplitudes. The dotted lines are the predictions 
of Dickinson et al. (2003) for the Ha-correlated free-free emission 
assuming values for the electron temperature between 5000 K 
and 9000 K (the lower and upper bounds respectively). There is 
weak evidence that the data prefer low values for f^, moreover if 
we consider our 'preferred' value for the electron temperature of 
7000 K it appears that is constrained to be less than ~ 0.35. 



then the 95% confidence limits on the scaling amplitude im- 
ply some joint constraint on fd and Te. By comparison of the 
predicted scaling of the Ha map, for the temperature range 
5000-9000 K (discussed in section|22l, with the 95% confi- 
dence range for the 31.5 GHz coupling parameters, we infer 
a weak limit on fd of ;6 0.75, but this also implies a value 
for Te lower than the 7000 K observed locally. Indeed, if we 
adopt this local value as the typical value for Te then fd is 
constrained to lower values still, < 0.35, in good agreement 
with Dickinson et al. (2003). 



4.4 Constraints from the 19.2 GHz data 

Here we apply the analysis of the previous sections to de- 
rive the correlation coefficients between the three foreground 
templates and the 19.2 GHz data of Cottingham (1987). 
This will allow us to place stronger constraints on the syn- 



chrotron and free-free emissions which rise with decreasing 
frequency and to investigate any dust-correlated component 
at lower frequencies. The results are summarised in Fig. El 
The analysis has again assumed a scale-invariant (n = 1) 
CMB component with normalisation Qrms-ps= 18 piK. We 
do not find any significant dependence on the choice of these 
parameters over the range allowed by the DMR data^. We 
also find that the spatially varying noise description is pre- 
ferred over the constant noise level per pixel used in the 
previous analysis by de Oliveira-Costa et al. (1998). 

Fig. El shows that there is a highly significant correla- 
tion between the 19.2 GHz data and all three foreground 
templates, particularly for the \h\ > 15° cut. The free-free 
coupling constant calculated as a free parameter agrees most 
closely with the expected emission at fd = 0; high values of 
fd are rejected. 

The stronger synchrotron signal at 19.2 GHz compared 
with that at DMR frequencies enables a more definitive spec- 
tral index between 408 MHz and these higher frequencies to 
be obtained as shown in Table 2. The data are consistent 
with a brightness temperature spectral index oi 13 = 3.1±0.1 
at these frequencies. 

The 140 /xm dust couplings are of lower amplitude than 
those determined in de Oliveira-Costa et al. (1998); their 
analysis did not have the Ha template available to it, and 
moreover they removed the observed quadrupolar contribu- 
tion to the sky signal before the correlation study. We also 
repeated our analysis in a similar fashion, using only the 
DIRBE 140 lira, and 408 MHz sky templates, and we indeed 
find results more consistent with this earlier treatment, al- 
though approximately Icr lower in amplitude. We also note 
that de Oliveira-Costa et al. (1998) found that their statisti- 
cal error was significantly lower than the dispersion in cou- 
pling amplitudes derived from correlations of the data with 
random rotations of the Galactic templates. The statistical 
error determined for our analysis is completely consistent 
with the value determined by such a sky rotation approach. 
We conclude that our method is unbiased, and affords sta- 
tistically robust error estimates. 

In the remainder of the paper, we will consider only 
those results where we have subtracted the predicted free- 
free contribution from the microwave data, and we adopt 
fd ~ and Te = 7000 K, the values which are maximally 
consistent with the picture we have established at this point. 



5 THE NATURE OF THE ANOMALOUS 
DUST-CORRELATED EMISSION 

We have shown that there is a diffuse foreground component 
in the COBE-DMR data at 31.5 and 53 GHz and the 19.2 
GHz data which is strongly correlated with the 140 fj,m dust 
emission. In order to make further progress in determining 
the nature of the anomalous component of dust-correlated 
emission, it is necessary to identify its spectral behaviour 
over the range of frequencies studied. This in turn requires 
that we separate out from the higher frequency observations 
that fraction of the observed dust correlation which is asso- 
ciated with thermal dust emission. 

This is consistent with the analysis of Boughn et al. (1992) 
which placed only upper limits on the normalisation amplitude. 
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Figure 3. Derived coupling coefficients between the 19.2 GHz 
data and three templates for foreground emission. Nomenclature 
and symbols as in Fig. 



There are several possible approaches to modelling the 
thermal dust at DMR frequencies. The first is to use an 
IRAS or DIRBE template of FIR emission and extrapolate 
to microwave frequencies using a simple model of the dust 
emissivity and temperature. In particular, we could use the 
model from DL98 where the emissivity index is set to ~ 1.7 
and the temperature and optical depth are adjusted so that 
the peak dust emission occurs at 140 ^m, or set the emissiv- 
ity index to ~2 and dust temperature to 17.5 K, as derived 
in Lagache et al. (2000). Kogut et al. (1996a) combined the 
observed correlations between the DMR and DIRBE 140 
^m data with additional information about r.m.s fluctua- 
tions obtained directly from the DIRBE 100, 140 and 240 
/im sky maps to attempt to infer constraints on the thermal 
dust emission. These results are quite consistent with the 
two previous models. Another approach, which we adopt 
here, is to use a model which incorporates data available 
over a wide range of FIR and microwave frequencies; suit- 
able models have been generated by FDS as described in the 
next section. 



5.1 Estimation and removal of thermal dust in 
the DMR data 

FDS have used 7° resolution data in 123 frequency lines be- 
tween 100 and 2100 GHz (A = 3 mm to 0.14 mm) from 
COBE-FIRAS data along with COBS-DIRBE data in the 
range 100 to 240 fj,m to produce models of dust emission 
which can be extrapolated to DMR frequencies. Their mod- 
els 7 and 8 afford the most acceptable fits to the observa- 
tions, according to a chi-square criterion. We have adopted 
model 7, differing little from model 8, which has two dust 
components, one at 9.6 K with ct = 1.5 and the other at 
16.4 K with a = 2.6. We proceed by subtracting from the 
DMR data the model 7 predictions at each frequency, then 
correlate the residuals against the DIRBE and Haslam tem- 
plates as before, with either an unconstrained and fitted free- 
free contribution, or after subtracting the predicted free- free 
component assuming fd~0 and Te= 7000 K. Fig.^Jshows 
the results from this analysis. The conclusions reached based 
on model 8 are indistinguishable from those presented here 
(although model 7 has ~ 20% higher r.m.s fluctuations at 
90 GHz relative to model 8). 

It is apparent when comparing Fig. 2] with the left- 
hand panel of Fig. that the thermal dust emission ac- 
counts for essentially the entire observed correlated sig- 
nal at 90 GHz. The best estimate of the dust components 
comes from the |6| > 15° cut and fixed Ha data. Here we 
find that at 90 GHz the model thermal dust component 
is equivalent to 1.39 /iK/(MJy sr~^) at 140 fim leaving 
0.77 ± 0.37 /iK/(MJy sr~^) of anomalous dust-correlated 
emission. This signal may be indicative of a small anoma- 
lous component, but equally it could be a consequence of the 
FDS parameterisation inadequately predicting an enhanced 
(very) cold dust component (for example, see Reach et al. 
1995). However, with the current data it is premature to 
radically alter the FDS model. 

The most important result of this study is the high 
statistical significance of the identified anomalous dust- 
correlated component of 6.06 ± 0.60 pK/MJy sr~^ at 31.5 
GHz. At 53 GHz the dust-correlated DMR coefficient is 
1.85±0.36 /iK/(MJy sr~^), of which the thermal dust model 
7 contributes 0.57 /xK/(MJy sr~^). The two dust-correlated 
components will therefore be approximately equal in ampli- 
tude at frequencies of ~ 60 GHz and suggests that this may 
be close to the minimum in terms of Galactic foreground 
contamination on large (~ 7°) angular scales. The dust con- 
tribution is discussed further in sections 15.41 & |7| 

5.2 Characterising the anomalous dust-correlated 
spectrum 

In order to quantify the nature of the anomalous dust- 
correlated component, we first compute the effective power- 
law spectral index between pairs of frequencies, as tabulated 
in TableEl We use the 19.2 GHz resuhs directly from Fig.gl 
without further correction for the thermal dust contribution 
since this is negligible at this frequency (comprising ~ 50% 
of the thermal contribution at 31.5 GHz which is in itself 
negligible). In addition, all data points are fitted simultane- 
ously by a model of the form Aanom (^) where Aanom, 
in units of /iK/(MJy sr~^) referred to the 140 /im template, 
is the amplitude of the emission at the reference frequency 
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Figure 4. Derived coupling coefficients between the DMR data at 
31.5, 53 and 90 GHz and three templates for foreground emission. 
Nomenclature etc as in Fig. Here, we remove the contribution 
from thermal dust emission by correcting the DMR, data before 
fitting the templates using the predictions and templates from 
FDS model 7. The free-free contribution is estimated only for the 
template derived assuming = 0. The results show that the 
majority of the emission at 90 GHz is due to thermal dust, with 
anomalous dust-correlated emission dominating at 31.5 GHz, and 
a mix of the two at 53 GHz. 



vo, which we take as 31.5 GHz. The values are consistent 
with a spectrum steeper than the typical value 2.15 of free- 
free emission, although in themselves not ruling out such an 
emission mechanism. Previous interpretations had proposed 
free- free as the source of the excess dust-correlated emission, 
yet in this work we have utilised a sky map of Ha optical 
line emission as an explicit tracer of this foreground com- 
ponent. The derived correlation with the data is consistent 
with that predicted for free- free emission and moreover, does 
not ehminate the dust-correlated emission. 

Microwave data on the anomalous dust-correlated emis- 
sion in addition to that given in the DMR and the 19.2 GHz 
surveys made at intermediate Galactic latitude {\b\ ^ 20°), 
are presented in Fig. |S] To allow a direct comparison with 
these previous studies, we have recomputed the DMR and 
19.2 GHz correlation coefficients at 100 fim. These are a 
factor of ~ 2 higher than the corresponding values for the 
140 /im correlations, reflecting the weaker dust emission at 
100 fim. Both the total dust (filled circles) and anomalous 



Table 3. The effective spectral index between pairs of frequencies 
for the anomalous dust-correlated component (after correction for 
free-free emission assuming /d = and for thermal dust emission 
using model 7 of FDS). The bottom line represents the best fit 

spectral index for a power law model of form Aanom (7^) 
where uq is taken as 31.5 GHz computed over all frequencies. 
68% c.l. errors are reported. 
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dust-correlated (open circles) emission are shown for com- 
parison. The best-fit power law is shown over the frequency 
range 19 — 100 GHz. Tenerife data on a 5°scale at 10 and 15 
GHz give correlation coefiicients of 71. 4± 14.2 and 51.4±8.0 
^K/(MJy sr"^) for \b\ > 20° at 100 nm respectively (de 
Oliveira-Costa et al. 2002). Mukherjee et al. (2001) sug- 
gested lower values based on a smaller Tenerife database. 
Finkbeiner et al. (2002) have also studied two isolated dust 
clouds on the Galactic plane at 5, 8.25 and 9.75 GHz; the 
dust emissivity differs by more than a factor of 10 in the two 
clouds and hence they have not been plotted. Nevertheless 
they find a dust-correlated emissivity which falls between 
10 and 5 GHz. There is a suggestion of a turn-over in the 
anomalous component spectrum at ~ 10 GHz but the ex- 
act shape is still poorly defined. Information at frequencies 
of 5 — 15 GHz is clearly important in defining the anoma- 
lous dust-correlated spectrum. A large area survey at in- 
termediate and high Galactic latitudes has been made with 
the Jodrell Bank Interferometer at angular scales of l°and 
2° (Jones et al. 2001). At 2°resolution the observed r.m.s 
level is 86 ± 12 ^K; this includes all the foreground compo- 
nents and the CMB. The synchrotron and free-free compo- 
nents comprise more than half this value, leaving ;S 50 fiK 
for dust emission. The corresponding dust emissivity at 5 
GHz is < 100 ^K/(MJy sr"^) at 140 fim. 

Fig. Elalso shows dust data at 150, 240 and 410 GHz 
from the BOOMERanG experiment (Masi et al. 2002) also 
for l&l > 20°, corrected from thermodynamic units to bright- 
ness temperature units. The FDS thermal dust model 7 nor- 
malised to 100 fim (dashed line) is in relatively good agree- 
ment with these points data. 

Fig. |S] therefore presents our assessment of the relative 
contribution of the anomalous and thermal components of 
dust emission in the range 10 to 500 GHz. The contributions 
are equal at ~ 60 GHz. These results apply to large areas 
of sky (iJ 1 sr) at intermediate and high latitudes where 
deep CMB observations are made; the situation closer to the 
plane near sources of heating and star-formation is likely to 
be different. 

The best-fitting model to the anomalous dust-correlated 
component is a power law model spectrum of the form f"'' 
with l3 ~ 2.5, in the range 19 — 90 GHz. This implies 
an emission component with a dust-like morphology but a 
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Figure 5. Dust emissivity correlation coefficients referred to 100 
/im, in units of /.iK/(MJy sr"'^) for |b| ^ 20°. Data points are: 
Tenerife at 10 and 15 GHz (filled triangles; de Oliveira-Costa et 
al. 2002), 19.2 GHz (filled square; this work), DMR at 31.5, 53 
and 90 GHz (filled circles; this work) and BOOMERanG at 150, 
240, and 410 GHz (filled stars; Masi et al. 2001). Also plotted 
are the DMR correlations after the thermal dust component has 
been subtracted (open circles) using model FDS7 (dashed line). 
The best fitting power-law model for the anomalous component 
is shown as a heavy dotted line over the range 19 — 90 GHz. 
Spinning dust models from DL98 for dust in 3 phases of the in- 
terstellar medium (CNM, WNM, WIM) are shown as dot-dashed 
lines scaled to 100 /im (see text). The preferred total spinning 
dust model from DL98 is shown as a solid line. 

synchrotron-like spectrum. The physical explanation of such 
a component is open to debate, and we explore dust-related 
alternatives below. It should be noted, however, that the ev- 
idence presented here for dust-correlated emission does not 
necessarily imply that the emission is intrinsic to the dust 
grains themselves. For example, it may arise from a source 
which is associated with the dust, such as may result from 
the star-formation process in a spiral arm environment. This 
is the interpretation taken by the WMAP team (Bennett 
et al. 2003) and is discussed in section ICTl 

5.3 Possible models of the anomalous emission 

A widely discussed model is in terms of rotational emis- 
sion from very small dust grains - a 'spinning' dust compo- 
nent, as proposed in DL98. They propose a three-component 
model of the spinning dust containing contributions from 
dust found in the Cold Neutral Medium (CNM), Warm 
Neutral Medium (WNM) or Warm Ionised Medium (WIM). 
The spectra for these three environments with the inferred 
anomalous emission are shown in Fig.|K|scaled to 100 /xm as 
dot-dashed lines. The preferred DL98 model has these com- 
ponents in the ratio 0.43, 0.43 and 0.14 respectively. The 
typical spectral index of the calculated dust emission over 
the frequency range 19 — 53 GHz is between 3.3 and 4, rea- 
sonably consistent with the observed values (Table|^. How- 
ever, the results presented in this paper indicate that, whilst 
the preferred spinning dust model of DL98 has a spectral 
shape which can reasonably match the data, some adjust- 
ment of amplitude is required. Fig.|K|shows the model (solid 



line) scaled to approximately fit the plotted data points. 
A more formal analysis indicates that an improved fit 
results from models in which emission from the WNM is al- 
most completely dominant. This can be understood rather 
easily: the CNM-component can be adjusted to give an ex- 
cellent fit to the data points at 19.2 and 53 GHz, but then 
over-predicts emission at 31.5 GHz; the WNM-component is 
scalable to provide good agreement with the data at 19.2 and 
31.5 GHz although then under-predicts the emission at 53 
GHz; the WIM-component can only be matched well to one 
frequency. Since the two lowest frequency data points have 
the most statistical weight, models containing a substantial 
contribution from the WNM are favoured. Spinning dust 
models therefore remain a possibility to explain the anoma- 
lous dust-correlated component, although the observed spec- 
trum is not as strongly peaked as the current dust models 
which therefore require some adjustments to match more 
closely with the new constraints. Given the uncertainties in 
the model predictions, depending on the composition, size 
and shape of the dust grains and conditions in the interstel- 
lar medium, this is certainly feasible. 

An alternative model for the anomalous dust-correlated 
emission at frequencies i, 100 GHz is due to thermal fluc- 
tuations in the magnetization of interstellar grains (Draine 
& Lazarian 1999) which results in magnetic dipole radia- 
tion. A population of Fe grains or grains with Fe inclusions 
are predicted to produce strong emission with a broad peak 
near 70 GHz due to the Frohlic resonance. Normal param- 
agnetic grains give much weaker emission at a level of 10~^ 
of the dipole emission expected from spinning dust grains. 
Hypothetical Fe-rich materials could produce radiation com- 
parable to that from spinning grains in the 14 to 19 GHz 
range. Draine & Lazarian (1999) conclude that magnetic 
dipole emission at v < 100 GHz might contribute significant 
emission if appropriate magnetic materials existed; both as- 
tronomical and laboratory measurements could clarify this 
situation. In the current context, the spectrum of the radi- 
ation is quite different; the magnetic dipole emission is rel- 
atively flat while the dipole (spinning) emission is peaked. 

A possible explanation for excess microwave emission 
in the North Celestial Pole region, at 14.5 and 31.7 GHz, 
was proposed by Leitch et al. (1997) in terms of free- free 
emission from ionized gas at Te > K. This explanation 
was considered unlikely because of the energetics in creat- 
ing this hot gas. Draine & Lazarian (1998a) argue that the 
case against this hypothesis is even stronger for the larger 
DMR sky coverage since the power radiated is a factor of 
~ 100 greater than that expected from the typical super- 
nova rate in the Galaxy. The current analysis also argues 
against this in view of the strong dust correlation and the 
fact that it persists in the presence of a simultaneous flt to 
the Ha template. However, it should be noted that the ratio 
of Hof to free-free radio continuum drops as the gas temper- 
ature is increased above ~ 3 x 10* K, hence the requirement 
for temperatures above ^ 10^ K. The deflnitive test of this 
model would be to analyse the correlation of microwave data 
with X-ray data at energies of ~ 1 keV. Indeed Kneissl et al. 
(1997) did find a significant (> 2a) correlation of DMR data 
with ROSAT all-sky X-ray data but only for the quadrupole 
£ = 2 component. After removal of the quadrupole term, no 
significant correlation remained. 

Another possible source of microwave emission from 
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dust grains is in terms of low-lying 2-level tunnelling tran- 
sitions at low temperatures. The laboratory measurements 
of Agladze et al. (1996) show surprisingly large and vari- 
able values for the millimetre (80 to 400 GHz) absorption 
coefficients of small silicate particles similar to interstellar 
grains; these effects occur at temperatures found in the in- 
terstellar medium (T J; 20 K). This phenomenon is also 
found in amorphous silica-based glasses and is believed to be 
due to resonant tunnelling absorption, although the detailed 
physics is not fully understood (Bosch 1978, Agladze & Siev- 
ers 1998). The low frequency spectrum of this emission has 
not yet been determined so a comparison with Fig. |^ is not 
possible. Further progress awaits laboratory measurements. 

5.4 Comparison with WMAP 

Since our more detailed re-analysis of the COBE-DMR 
data, the first results from the analysis of one year of ob- 
servations by the WMAP satellite have been released (Ben- 
nett et al. 2003a and references therein). Their foregrounds 
analysis (Bennett et al. 2003b) uses a maximum entropy 
method applied to the 5 frequency channels of the WMAP 
observations, together with data from outside the microwave 
spectral regime, specifically the Haslam 408 MHz survey, 
the full-sky Ha map published by Finkbeiner (2003) and 
a dust template derived from COBE-DIRBE and FIRAS 
data (FDS model 8 normalised at 94 GHz). These exter- 
nal data-sets are claimed to be the key to separating the 
three emission mechanisms (free-free, synchrotron and spin- 
ning dust) with a spectral index of /3 ~ 2 at the WMAP 
frequencies. One should note, however, that no unique spin- 
ning dust template has been employed since this putative 
emission falls off rapidly outside the microwave regime, and 
therefore there is no obvious candidate (although DL98 have 
proposed that the DIRBE 12 /im map might be used). More 
importantly, the analysis is allowed only to produce a com- 
bined synchrotron/spinning dust solution at each frequency, 
with no attempt made to separately disentangle these two 
components. As a result of the analysis, they identify a dust- 
correlated component in the lower frequency (23, 33 and 41 
GHz) WMAP channels. However, they interpret this as a 
synchrotron component with a spectral index of /3 ~ 2.5, a 
value which we emphasise is in close agreement to that de- 
rived in this analysis using COBE-DMK data. They argue 
that this emission originates in star-forming regions close 
to the plane while the steeper (/3 ~ 3) synchrotron compo- 
nent is from gas in the halo, a picture which is consistent 
with what has been found in the integrated radio and FIR 
emission for several nearby star-forming galaxies. 

In their figure 9, Bennett et al. (2003b) show the steep- 
ening of the mode for the spectrum of the synchrotron plus 
spinning dust component at 23-41 GHz and conclude that 
this is evidence for a spectral break, corresponding to a 
steepening of A/3 ~ 0.5 at frequencies of ~ 20 GHz. How- 
ever, this could also be the signature of spinning dust. Cur- 
rent models of spinning dust (DL98) also give a steepening 
of index with spectral indices ranging from ~ 2.5 at 30 GHz 
and steepening to greater than 4 at frequencies 50 GHz. 

An independent re-analysis of the WMAP foregrounds 
has been performed by Lagache (2003). WMAP data is cor- 
related with the neutral Hi component and an excess is simi- 
larly found at longer wavelengths. By comparing the correla- 



tions with the Hi column density, there is evidence that the 
excess emission is associated with small transiently heated 
dust particles. The exact physical mechanism of such emis- 
sion remains to be determined, but possible models certainly 
include the spinning dust rotational emission mechanism or 
Very Small Grain (VSG) long-wavelength emission resulting 
from transient heating by UV photons. 

In the previous discussion, we have attempted to 
present arguments as to why the interpretation of the 
WMAP foreground results remains open, and that a spin- 
ning dust component can be compatible with the observa- 
tions. However, we would like to emphasise that we do not 
consider that this is an unambiguous solution to the origin 
of the dust-correlated emission. The results of our analysis 
based on independent data sets are completely consistent 
with that of the WMAP team, and indeed the flat spec- 
trum fit we have derived to the anomalous dust-correlated 
component is formally better than any spinning dust model 
considered. In order to clearly distinguish the various mod- 
els for this excess emission, low frequency data (few GHz 
to 15 GHz) will be vital. Evidence for WMj4P-consistent 
spectral indices determined from lower frequency radio mea- 
surements may be important to corroborate the star-forming 
region, flat spectrum hypothesis. We await the forthcoming 
publication of a full-sky 1420 MHz survey with interest. Fi- 
nally, detailed spatial comparisons between the radio and 
FIR emission will be critical in determining the physical 
processes associated with this dust-correlated component. 



6 RESULTS - CMB 

We are now in a position to consider the implications of 
our better understanding of the Galactic diffuse foregrounds 
upon estimates of Qrms-ps, the CMB quadrupole normal- 
isation and n, the power-law spectral index. In section |3 
we describe how the simultaneous solutions for the cou- 
pling coefBcients of the DMR data to the three foregrounds 
forces the CMB distributed power to be invariant between 
the three DMR frequencies, so that each set of foreground 
solutions has an accompanying estimate of Qrma-ps and n. 
Table 0| gives the Qrms-PS and n values for the following 
areas of parameter space. 

• Column 1: Foregrounds used in the fit- 140 /im dust 
(D), synchrotron (S), free- free (F), free- free forced to the 
Ha-derived value (FixedF) and a flxed thermal dust compo- 
nent following FDS model 7 (FDS7). 

• Column 2. The Galactic cuts used in the flts- straight 
Galactic latitude cuts at |fo| > 15° and 30°, and the standard 
DMR cut at |6| > 20° which also excludes extra regions of 
stronger emission. 

• fd, the correction applied to the free-free (Ha) template 
for foreground dust absorption. We consider fd = 0.0 is the 
best value for the latitudes considered here (see section 

An examination of Table 2] immediately shows that the 
solution for Qrms-PS and n is strongly independent of the 
parameters incorporated in the present study. We believe 
that our best estimate of Qrms-PS and n are given for 
fd = 0, fixed free-free and the FDS dust model included, 
namely Qrma-ps= 15.18l2'3o ^K, n = 1.231q'23- This may 
be compared with Qrms-ps= 16.4912 49 A'K, n = 1.19_q'23 
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if only the dust and synchrotron foregrounds are used with 
an aggressive Galactic cut; the earlier analysis of Bennett et 
al. (1996) gave Qrms-ps= IS.Sla g ^K, n = 1.2 ± 0.3. The 
present analysis gives values in excellent agreement with the 
earlier values, but with smaller errors. 

The reason that the solutions for Qrms-ps and n in 
Table are not strongly dependent upon the precise values 
of the Galactic foregrounds used is that the weighting of the 
data from the three radiometers with respect to the likeli- 
hood analysis is strongest for 53 and 90 GHz; the weightings 
are approximately 8, 65, and 27 per cent at 31.5, 53 and 90 
GHz respectively, which reflects the corresponding noise lev- 
els of the three frequencies. As will be shown in section |7| 
the minimal foreground contributions are in the range 50 to 
100 GHz and therefore explains why there is only a small 
effect on Qrms-PS and n. 

Indeed, even for a Galactic cut of 15°, there appears to 
be remarkable consistent with previous analyses, provided 
the three templates for foreground emission are included 
in the analysis. The worst perturbation in Qrms-PS and 
n corresponds to the result from the data corrected for dust 
and synchrotron alone for this cut. However, even here the 
changes are less than la for the quadrupole normalisation, 
and a less significant amount for the cosmological spectral 
index value. Clearly, the uncorrected fit for this latitude cut 
is inconsistent at worse than 2a. 

Given the success of the Galactic correction, we might 
wonder whether it is possible to be yet more aggressive, 
and attempt to correct for Galactic contamination at still 
lower latitudes. We have tried this using a cut in which we 
only eliminate those pixels lying within 10° of the Galactic 
plane, but in this case we find that, whilst the derived cor- 
relations of the foreground templates with the DMR maps 
are very consistent with the 15°results discussed above, the 
likelihood analysis begins to deteriorate. This close to the 
Galactic plane, there remains a sufficiently complex struc- 
ture in the Galactic emission that it is inevitable that simple 
correlation studies which do not allow for variations in spec- 
tral behaviour cannot succeed in subtracting the foreground 
at the necessary level of accuracy. Moreover, for this case, 
half of the 7° FWHM beam area of DMR lies at |6| < 10°. 
Residual poorly corrected features contribute to an overall 
flattening of the cosmological power-law spectral index and 
increase in normalisation amplitude well beyond the allowed 
ranges yielded by the studies confined to more conventional 
latitude cuts. It is unlikely that even sophisticated analy- 
sis methods, such as those based on independent compo- 
nent analysis (Maino et al. 2002) or Maximum Entropy in 
harmonic space (Stolyarov et al. 2002) will be able to cir- 
cumvent these issues. With higher sensitivity, higher angular 
resolution observations such as those due to WMAP and ex- 
pected for Planck, such methods as applied to small regions 
or even on a pixel-by-pixel basis may yield an improved un- 
derstanding of the physical nature of the foreground emis- 
sion and allow greater sky coverage to be used for cosmo- 
logical studies. 

7 CONCLUSIONS 

In this paper we have reappraised the contribution of Galac- 
tic foreground emission to the COBE-DMK data, aug- 



mented with lower frequency information at 19.2 GHz. It 
has been demonstrated that a new template for Galactic 
free-free emission derived from observations of the Ha emis- 
sion is in excellent agreement with expectations, and further- 
more allows us to claim, with reasonable certainty, that the 
amount of emission absorbed by foreground dust is low at 
high Galactic latitudes, in particular for the preferred local 
value of the electron temperature of 7000 K. After subtrac- 
tion of the free-free contribution, there is continued and un- 
ambiguous evidence of anomalous dust-correlated emission 
at low frequencies, consistent with previous analyses. Sub- 
traction of an FDS model for Galactic thermal dust emission 
indicates that the correlated dust component at 90 GHz is 
essentially dominated by a vibrational (thermal) emission 
mechanism. The anomalous component at lower frequencies 
(19 — 90 GHz) has a power-law spectral index of ~ 2.5, 
somewhat flatter than the putative spinning dust compo- 
nent proposed in DL98. Indeed, their preferred mixture of 
spinning dust emission from a three-component dust model 
requires some adaptation to afford an improved fit. The level 
of Galactic synchrotron emission is consistent with expecta- 
tions based on lower frequency observations at 408 MHz and 
a steep temperature spectral index of ~ 3.1. The best-flt cou- 
pling coefficients (|fe| > 15°) for each foreground component 
are given in Table for 19.2, 31.5, 53 and 90 GHz. Both 
the total and anomalous dust-correlated are given for com- 
parison while the free-free scalings are those predicted by 
the Ha template assuming Te = 7000 K and fd ~ 0. These 
results are completely consistent with the recent WMAP 
analysis and independently confirm their results on large 
angular scales. In particular, we demonstrate the existence 
of an emission mechanism that has a spatial morphology 
similar to dust but with a synchrotron-like spectrum in the 
range 19-90 GHz. 

Our better understanding of the foregrounds allows us 
to correct for them over a larger region of the sky and re- 
assess the value of Qrms-PS and n. The most signiflcant 
changes in general result from an analysis made with an 
aggressive |fe| > 15° Galactic cut, yet even here the results 
are compatible with earlier results at better than the ~ la 
level, although the error is slightly smaller, suggesting some 
possibility of making use of modestly larger sky coverage for 
the inference of cosmological parameters. Interestingly, what 
might be considered our best model for the foreground cor- 
rection - with this aggressive Galactic cut, utilising the FDS 
Model 7 dust template to correct for thermal dust emission, 
the Ha template from Dickinson et al. (2003) to trace the 
free- free emission (flxed assuming that the electron tempera- 
ture is 7000 K), together with unrestricted fits to the DIRBE 
140 fim and 408 MHz sky surveys to correct for anomalous 
dust-correlated and synchrotron emission respectively - re- 
sults in shifts in Qrms-PS and n of only ~ 0.6 /iK and ~ 0.01 
relative to the previous standard uncorrected result. All re- 
sults with fd = are consistent with values of Qrms-PS ~ 
15.3 and n — 1.2. 

Fig. m summarises the integrated amplitudes of Galac- 
tic foregrounds at large (^ 7°) angular scales over the most 
important frequencies for CMB studies. Our best models for 
synchrotron {l3 = 3.02 scaled to the 408 MHz map), free- free 
(scaled Ha predictions assuming Te = 7000 K; dot-dashed 
line) and thermal dust model (FDS7; dashed line) emissions 
are plotted in terms of thermodynamic temperature units 
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Table 4. Cosmological parameter fits to the DMR data for a power law model parameterised by r.m.s quadrupole normalisation amplitude 
Qrms—PS (mK) and power law index n. The fits are made for three Galactic cuts, and after correction for various combinations of 
foreground emission as determined by simultaneous fits to template sky maps. The corrections should be understood as follows: D - 
dust traced by the DIRBE 140 fim sky map; F - free-free traced by the Ha template of Dickinson et al. (2003) after correction for dust 
absorption specified by = 0, 0.5 and 1.0; S - synchrotron emission determined from the 408 MHz survey of Haslam et al. (1982); 
FixedF - we subtract a fixed free-free contribution determined by scaling the Ha template using the amplitudes predicted in Dickinson 
et al. for an electron temperature of 7000 K after correcting for dust absorption specified by = 0, 0.5 and 1.0; FDS7 - we subtract a 
thermal dust contribution as predicted by FDS Model 7. 



Correction 


Coverage 




Qrms-PS 






n 








fd = 0.0 


0.5 


1.0 


fd =0.0 


0.5 


1.0 


Uncorrected 


b > 15° 
standard 
b > 30° 


24.45ll«° 
15.85l^:fg 
15.3ll^;|5 


- 


- 


Q 94+0.21 
"-"^-0.22 

-I 99+0.24 

1-^^-0.26 

1 20+°-28 
1-^^-0.29 


- 


- 


DS 


b > 15° 
standard 
b > 30° 


16.491^;!^ 
15.081^:°^ 

15.02t3;39 






1 in+0.22 
1-1^-0.23 
-I 90+0.24 
l-^'-'-0.25 
1 ig+0.27 
l-l='-0.29 






DSF 


b > 15° 
standard 
b > 30° 


15 34+2-91 

14 70+3-39 
It. (U_2 g2 


15.1812-84 

15 07+3.10 

iO.O(_2 49 
15.18_2;72 


1 5 1 2+2-81 
1^-1^-2.30 
15 34+3.10 
10.0t_2 4g 

15-18-2^72 


-I 99+0.22 
1-^^-0.24 
-I 91+0.24 
1-^1-0.25 
1 99+0.27 
1-^^-0.29 


1 90+0.22 
l-^'^-0.23 
1 20+°-24 
1 10+0. 27 
l-l''-0.29 


-I 90+0.22 
l-^'^-0.23 

1 20+°-24 
1-^^-0.25 
-1 10+0. 27 
l-l*'-0.29 


DS-FixedF 


b > 15° 
standard 
b > 30° 


15.34l2.9J 

15.47l^:« 
15.6613-52 


15 24+2-91 
iO.Z4_2 30 

15.63+11 
15.56_2 7g 


15.341-931 

15.79l3;20 
15.6613:52 


1 99+0.22 
1-^^-0.24 

1 20+0-24 
l-^l'-0.25 
1 15+0-27 
l-l'^-0.28 


1 99+0.22 
1-^^-0.24 
1 10+0. 24 
l-l*-0.25 
1 15+0.27 
l-l''-0.29 


1 90+0.22 
l-^'^-0.24 
1 17+0.24 
l-l'-0.25 
1 15+0-27 
l-l"'-0.29 


DS-FDS7 


b > 15° 
standard 
b > 30° 


15 12+2-^4 
10.1Z_2 30 

15 18+3-07 

14 fi4+3-36 
lt.D'i_2 g2 


14.8612-75 

15.18l|«^ 
15 05+3-*2 

iO.UO_2 gg 


14.80l2;^« 

15-181^4^ 
15.05l3:« 


1 90+0.22 
l-^'^-0.23 
1 91+0.24 
1-^1-0.25 
1 99+0.27 
1-^^-0.29 


1 9cr+0.22 
l-^°-0.23 
1 91+0.24 
1-^1-0.25 
1 19+0-27 
1-1^-0.29 


1 26+°-22 
l-^'-'-0.23 
1 91+0.24 
1-^1-0.25 
1 19+0-27 
l-l''-0.29 


DS-FixedF-FDS7 


b > 15° 
standard 
b > 30° 


I5.l8l2;«3j 
15.28+11 

15.50_2 72 


14.86l2;75 
15.37l|i3 


14.83l2;8i 
15.5613:10 
15.5313:52 


1 90+0.22 
l-^'^-0.23 
1 20+0-24 
l-^l'-0.25 

1 16+0-27 
1-10-0.29 


1 9cr+0.22 
l-^''-0.23 
1 20+0-24 

1 16+0-27 
l-l"-0.29 


-, 97+0.22 
l-^'-0.24 
1 -10+0. 24 
1-1^-0.25 
1 15+0-27 
l-l"'-0.29 



Table 5. Definitive Galactic Template Coupling Coefficients; derived assuming a fixed free-free contribution with Te = 7000 K and 
fd = for the Galactic latitude range |6| > 15°. The FDS7 model for thermal dust emission has been applied in order to allow the 
anomalous dust-correlated emission contribution to be unambiguously determined. The Galactic coefficients have units /iK where 
X are the natural units of the template map: MJy sr~i for the 140 fivci DIRBE dust map, Kelvin for the 408 MHz synchrotron map 
and Rayleigh for the Ho (free-free) map. The free-free scaling coefficients are determined according to the relations in Dickinson et al. 
(2003), assuming = 7000 K for the central value, and 5000 and 9000 K for the lower and upper error bars respectively. 
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0.91 ± 0.83 
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-0.26 ± 0.56 
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2.16 ± 0.37 
0.77 ± 0.37 
-0.35 ± 0.54 
56+0-10 



(relative to Tcmb = 2.73 K). The best-fit anomalous compo- 
nent is shown (heavy dotted line) over the range 19—90 GHz. 
Data points are shown for synchrotron at 19.2 GHz (open 
triangle) anomalous dust at 19.2, 31.5, 53 and 90 GHz (open 
circles), total dust at 90 GHz (filled circle). The CMB sig- 
nal is plotted at a level of AT/T = 10"^ K. The approxi- 
mate total foreground level is represented by adding the four 
components in quadrature, with the anomalous component 
set to 52 fiK (the value at 19.2 GHz) for frequencies be- 
low 19.2 GHz. This suggests that the minimum foreground 



contamination for these data, is ~ 70 GHz. This is in good 
agreement with the WMAP analysis which predicts the min- 
imum to lie at ~ 60 GHz. 

The amplitudes for each of the foreground signals are 
now relatively well known as shown in Fig. |H] Strictly, Fig. 
|S| applies only for angular scales of 7° and larger, and for 
all-sky observations with a Galactic cut of jbj ^ 20°. For 
different Galactic cuts, there are likely to be variations of a 
factor of 2 — 4 for reasonable Galactic cuts as indicated in Ta- 
ble 1. Observations covering smaller regions of sky, perhaps 
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Figure 6. Foreground contaminations for the standard Galac- 
tic cut ^ 20°) at angular scales ^ 7°. Best- fitting models to 
19.2 GHz and DMR data are shown for each foreground compo- 
nent: synchrotron (dotted line), free-free (dot-dashed line), ther- 
mal dust (dashed line) and anomalous dust (heavy dotted line 
between 19 and 90 GHz). Units have been converted to ther- 
modynamic temperature fluctuations relative to Tcmb = 2.73 K; 
CMB fluctuations are represented as a horizontal line at a level of 
AT/TcMB = 10~^ for comparison. Data from this paper are: syn- 
chrotron 19.2 GHz (open triangle), anomalous emission at 31.5, 
53 and 90 GHz (open circles) and total dust emission at 90 GHz 
(filled circle). The thick solid line gives an estimate total fore- 
ground level by adding all the foreground models in quadrature. 
The frequency channels for the WMAP (short vertical dashed 
lines) and Planck (short vertical solid lines) GMB satellites are 
shown at the top of the figure. The Low and High Frequency In- 
struments (LFI and HFI) for Planck are drawn as heavy and fine 
lines respectively. 

cleaner regions (for example, away from the North Polar 
Spur), may be have reduced foreground contamination by 
a factor of ~ 2. Finally, at higher i"- values (smaller angular 
scales), it is known that diffuse Galactic foregrounds reduce 
in amplitude. The power spectra of Galactic emission have 
a typical power-law indices (Ci oc l'") of ~ 2 — 3 (Gia- 
rdino et al. 2001). Thus at 1° angular scales, the Galactic 
foregrounds, in terms of AT, will be reduced by a factor 
~ 2. However, at angular scales 1z 1°, discrete extragalactic 
sources become the dominant foreground. 

With relevance to the satellite missions WMAP and 
Planck, our analysis implies that the lowest level of fore- 
ground contamination from diffuse Galactic foregrounds at 
least on a 7° angular scale will be between 50 and 100 GHz. 
This has been confirmed by Bennett et al (2003b) and seems 
to hold to yet smaller scales ~ 1°. Our study also pre- 
dicts that for the WMAP escpeximent, the predominant fore- 
ground at 22, 30 and perhaps 40 GHz, will be due to the 
dust-correlated component. This is indeed found to be the 
case by the recent WMAP analysis, although the interpre- 
tation of the physical origin of the emission remains open 
to debate. More accurate determination of its spectral be- 
haviour, allied to lower frequency observations (5-15 GHz) 
and higher angular resolution observations, should allow the 
mechanism behind this unexpected foreground to be ascer- 
tained. Moreover, given the likely region of minimal fore- 
ground contamination in frequency space, the importance 



of the Planck-h¥\ instrument to the Planck mission as a 
whole cannot be overemphasised. 
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